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("^ , The ground state spin of the negatively charged nitrogen-vacancy center in diamond has been 

rS] ■ the platform for the recent rapid expansion of new frontiers in quantum metrology and solid state 

quantum information processing. In ambient conditions, the spin has been demonstrated to be a 
high precision magnetic and electric field sensor as well as a solid state qubit capable of coupling 
with nearby nuclear and electronic spins. However, in spite of its many outstanding demonstrations, 
r--s ' the theory of the spin has not yet been fully developed and there does not currently exist thorough 

*vj ' explanations for many of its properties, such as the anisotropy of the electron g-factor and the 

existence of Stark effects and strain splittings. In this work, the theory of the ground state spin is 
fully developed for the first time using the molecular orbital theory of the center in order to provide 
fj ■ detailed explanations for the spin's fine and hyperfine structures and its interactions with electric, 

C/3 , magnetic and strain fields. 

PACS numbers: 31.15.xh; 71.70.Ej; 76.30.Mi 
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I. INTRODUCTION 



The negatively charged nitrogen- vacancy (NV") cen- 
ter is a unique defect in diamond that has many promis- 
ing applications in quantum metrology and quantum in- 
formation processing (QIP). In particular, the ground 
state spin has been used in recent demonstrations of high 
precision magnetioi"— and electric field^ sensing, as well 
as spin-photon^ and spin-spinii"— entanglement. The 
NV~ center has also been employed to explore the devel- 
oping decoherence based sensing techniques'^"— Each of 
these demonstrations exploit the interaction of the spin 
with some configuration of electric, magnetic and strain 
fields and the center's remarkable capability of optical 
spin-polarization and readout i^i^ The demonstration of 
electric field sensing)^ which required the precise control 
of the spin using magnetic fields and the intricate mod- 
eling of the spin's interaction with electric, magnetic and 
strain fields, highlighted the NV~ center as a universal 
field sensor at the nanoscalc as well as emphasized the 
requirement for a more detailed theoretical understand- 
ing of this important spin in diamond. The development 
of such an understanding will enable more precise control 
and modeling of the spin in its current applications and 
also provide the necessary insight to explore its future 
applications. 

The NV~ center is a point defect of Cst, symmetry in 
diamond consisting of a substitutional nitrogen atom ad- 
jacent to a carbon vacancy (refer to Fig. [IJ. The center's 
electronic structure is summarized in Fig. [5] It consists 
of a '^^2 ground triplet state, an optical '^E excited triplet 
and several dark singlet statesi^ The fine structure of the 
^E excited triplet is highly dependent on temperature^^ 




FIG. 1: (color online) Schematics of the nitrogen-vacancy 
center and lattice depicting the vacancy (transparent), the 
nearest-neighbor carbon atoms to the vacancy (black), the 
substitutional nitrogen atom (brown), and the next-to-nearest 
carbon neighbors to the vacancy (white). The adopted co- 
ordinate system and carbon labels are depicted in the right 
schematic. 



and crystal strain)^ whereas the fine structure of the 
■^^2 ground triplet state is observed to be only weakly 
dependent on temperature^^ and crystal strain with a 
single zero-field splitting of Dgs ~ 2.87 GHz between the 
rris = and rus = ±1 spin sub-levels. At ambient tem- 
peratures, the fine structure of the excited triplet state 
replicates the ground triplet state with a single zero-field 
splitting of Des ^ 1-42 GHa^^i^ independent of crystal 
strain due to the dynamic Jahn- Teller effecti^ Zeeman 
and Stark splittings have been observed in the fine struc- 
tures of both triplet states,— i^Si"— although the Stark 
effect in the ground triplet state is several orders of mag- 
nitude smaller than in the excited triplet state i^^i'^° 
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FIG. 2: The electronic orbital structure (left) and fine struc- 
ture (right) at ambient temperatures formed from the ground 
afe^ and first excited aie^ molecular orbital configurations. 
The observed optical zero phonon line (1.945 eV)^ and in- 
frared zero phonon line (1.190 eV)^ transitions are depicted 
as solid arrows in the orbital structure. 



Electron spin resonance (ESR)^ and ab initio 
studies^i"— have confirmed that the electronic states 
of the center are highly localized to the vacancy and 
its nearest neighbors. The localization of the center's 
electronic states supports the application of a molecular 
model, in which the center's states arc approximated by 
configurations of molecular orbitals (MOs). The molec- 
ular model has been successfully applied to describe the 
effects of electric, magnetic and strain fields on the fine 
structure of the excited triplet state,— i^i^ however the 
model has not yet been applied to describe the effects of 
the fields on the ground triplet state. This has been due 
to an absence of the spin-orbit and spin-spin induced 
couplings of the center's electronic states, which have 
been obtained recently in Ref. [22|. Consequently, until 
now, the measurements of the effects of the fields on the 
ground state spin have been interpreted using the effec- 
tive spin-Hamiltonian formulism of ESRi^i^i^i^ Whilst 
the spin-Hamiltonian formulism has provided a practical 
model for the implementation of the center's applications 
to date, it docs not facilitate the correlation of the prop- 
erties of the ground state spin to the other properties of 
the center nor the theoretical prediction of more intricate 
properties of the spin. Consequently, the application of 
the molecular model to describe the effects of electric, 
magnetic and strain fields on the ground state spin is 
expected to provide the required detailed theoretical un- 
derstanding of the spin. 



In this article, the well established molecular 
model22ii^i^'^ of the NV~ center will be applied in or- 
der to fully develop the theory of the ground state spin. 
The fine and hyperfine structure and their corresponding 
eigenstates will be constructed prior to examining the ef- 
fect of electric, magnetic and strain fields on each. By 
using the matrix representations derived in the recently 
published electronic solution;^ explicit expressions in 
terms of the center's MOs will be derived in this work for 
the hyperfine interaction with the ^^N nucleus, the com- 
ponents of the electron g-factor tensor, and the Stark and 
strain interactions. The derivation of these expressions 
enables the rigorous definition of the spin-Hamiltonian 
of the ground state spin and the correlation of the ac- 
curately measured parameters of the spin to the other 
observed properties of the center. The expressions will 
also assist future ab initio studies to independently cal- 
culate the properties of the ground state spin. 



II. ELECTRONIC FINE STRUCTURE AND 

INTERACTIONS WITH ELECTRIC, MAGNETIC 

AND STRAIN FIELDS 

By adopting the adiabatic approximation and consid- 
ering the nuclei of the crystal to be fixed at their equilib- 
rium coordinates Rq corresponding to the ground elec- 
tronic state, the electronic Hamiltonian He of the NV^ 
center can be defined as 

He = J2^' + VNe{n,Ro) + Vso{x^,Ro) + f,./(x„ Xq) 

i 

+ ^ KeCXi, Xj) + Vss{Xi,:>ij) (1) 

where x.^ = (f^, Si) denotes the collective spatial and spin 
coordinates of the i*'' electron of the center, Xo = {Rq, I) 
denotes the collective equilibrium spatial and spin coor- 
dinates of the crystal nuclei, Ti is the kinetic energy of 
the i*'' electron, Vnc is the effective Coulomb potential of 
the interaction of the nuclei and lattice electrons with the 
electrons of the center, Vso is the electronic spin-orbit po- 
tential, Vhf is the hyperfine potential of the interactions 
between the crystal nuclei and the electrons of the center, 
Vee is the Coulomb repulsion potential of the electrons of 
the center and ^^s is the electronic spin-spin potential. 

Ab initio studies^i"— have confirmed the presence of 
three MOs (ai, Bx-, Cy) in the band gap of diamond 
and the center's observable electronic structure has been 
shown to consist of the ground a^e^ and first excited 
flie'^ MO configurations formed from the occupation of 
the MOs by four electronsi^ Note that the other two 
electrons of the six electrons associated with the center 
occupy delocalized Ai symmetric MOs within the dia- 
mond valence band and do not influence the observable 
properties of the center. The electronic states can be 
constructed by firstly defining orbital states with well de- 
fined Csi, orbital symmetry formed from products of the 



MOs and defining spin states with well defined Csy spin 
symmetryi^ Secondly, electronic states $^° j. that trans- 
form as specific rows k of irreducible representations j of 
the C31, group in spin-orbit space are formed from linear 
combinations of the orbital and spin state productsj^ 
Note that the quantum number n denotes the fine struc- 
ture level of the electronic state. For example, using the 
irreducible representations contained in Ref. 1471 the elec- 
tronic states of the ground triplet are 
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where $^2 = -75(01016x6^ — aiaieyCx) is the A2 orbital 
state of the ground triplet, the kets on the right hand 
side denote Slater determinants (overbar denoting down- 
spin), and the symmetrised S = 1 spin states in terms 
of the Sz eigenstates {jS*, m,)} arc 5^2 = |l:0)j '^e x = 
^(|1,1) + |1,-1)) and SE,y = ^(|l,i> - |1,-1». ' 

The majority of the spin-orbit states $^° ■ j. are eigen- 
states of the orbital components of the electronic Hamil- 
tonian Ho ^J2i'^t + VNein, Rq) + J2i>j Ke(xi, xj) with 
orbital energies denoted by Ej-s (where to remain con- 
sistent with Ref. [22, J denotes the irreducible represen- 
tation of the orbital state and 5* denotes the total spin 
of the spin states used to form $^° ■ j, and the orbital en- 
ergy of the ground triplet is defined as Ea2;1 = 0). Only 
the spin-orbit states ($3°b,x. ^3°B,r *i?£;.^' *9°B,y) asso- 
ciated with the ^E and ^E' singlets arc not eigenstates of 
Ho and are mixed by the Coulomb coupling coefficient n 
(refer to Ref. [22| for further details). The energies Ej-s 
including the effects of the Coulomb repulsion of the E 
singlets form the orbital structure of the center depicted 
in Fig. H 

The electronic spin-orbit and spin-spin potentials can 
be treated as first-order perturbations to Ho using the 
orbital energies Ej-s and the spin-orbit states $^° f, (ac- 
counting for Coulomb coupling) as the zero-order ener- 
gies and states of the perturbation expansion. The per- 
turbed energies correct to first-order En have been shown 
to be consistent with the observed fine structure of the 
center depicted in Fig. [5)^ It is found that the fine 
structure of the ground triplet state is governed by elec- 
tronic spin-spin interaction which splits the m^ — and 
nis = ±1 spin sub- levels such that £'2 — -Ei = Dg 
GHz, where 
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^0 is the vacuum permeability, ge = 2.0023 is the free 
electron g-factor, fiB is the Bohr magneton, fi = XiX + 
Tjiy + ZiZ {x,y,z being unit coordinate vectors), fi2 = 
^2 — n., and Z12 = Z2 — z\. As obtained in Ref. |22| . 



TABLE I: The spin coupling coefficients of tlie ground triplet 
correct to first-order in spin-orbit and spin-spin interactions. 
The Coulomb coupling coefficient k, the spin-orbit parameters 
(A||, Ax), the spin-spin parameters (-Di,£;,i,-Di,b,2), and the 
spin-spin coupling coefficient of the excited triplet t\ are as 
defined in Ref. [13 Note that iV^ = (1 
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the first-order corrected spin-orbit states $^°'- j. have the 
general form 
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first-order spin coupling coefficients. Using the results of 
Ref. [23, the first-order spin-orbit states of the ground 
triplet are 
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where k = x,y and the spin coupling coefficients of the 
ground triplet are contained in table |T] and are functions 
of the orbital energies Ej.g, the axial Ay = 5.3 GHz^ 
and non- axial Aj^ ^ GHz spin-orbit parameters, the spin- 
spin parameters Diea '^ MHz and Di^e,2 '^ MHz, the 
Coulomb coupling coefficient k, and the spin-spin cou- 
pling coefficient of the excited triplet 7/ = 0.053.— 

Note that recent strain measurements of the in- 
frared zero phonon line (ZPL)^ have indicated that the 
Coulomb coupling coefficient is significant k ~ 0.3 and 
therefore must be retained to second-order in the spin 
coupling coefficients. Furthermore, it should be noted 
that since the spin-orbit and spin-spin parameters are 
expected to be of the order of --10"^-10~^ eV (MHz- 
GHz) and the orbital energies are expected to be of the 
order of ~ 10~^ — 10"'^ eV, the spin coupling coefficients 
are expected to be of the order of 10~^ — 10^^*^. The rel- 
ative magnitudes of the different coefficients will become 
important in determining the leading order terms that 
contribute to the interactions of the ground state spin 
with electric, magnetic and strain fields. 

Given ([5]) and the zero-order orbital and spin operator 
matrix representations contained in Ref. l22l matrix rep- 
resentations of the ground triplet using the basis of first- 
order corrected spin-orbit states {^{°'a^,^2°'e x^'^T'e y} 



can be constructed for general one-electron orbital ten- 
sor operators Opg = J^t Op^qiji), which transform as the 
row q of the irreducible representation p of the C^v group, 
and the total spin operator S = X^i Si (see tabic |n|). The 
matrix representations are expressed in their most sim- 
plified form in terms of one-electron reduced matrix el- 
ements and the center's MOs. The one-electron matrix 
elements and the associated reduced matrix elements are 
related by the Wigner-Eckart theoren?^ 
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where cfij^k and 0/_g are MOs of symmetry (7, k) and (/, g) 
respectively, and (j) are the Clebsch-Gordon coefficients 
defined in Ref. [43. The detailed model of the inter- 
actions of the ground state spin with electric, magnetic 
and strain fields can thus be developed by applying the 
matrix representations of tabic |lT] to each interaction in 
turn. As a result, the interactions will be expressed in 
their most simplified form in terms of the spin coupling 
coefficients, the MOs, and the reduced matrix elements. 



A. Interactions with magnetic fields 

Defining B to be the applied magnetic field that is 
assumed to be approximately constant over the dimen- 
sions of the NV~ center, the interaction of the center's 
electrons with the magnetic field is described by the 
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where I = IxX + lyy + IzZ = r x p is the electron orbital 
magnetic moment, pis the electron momentum, h is the 
reduced Planck constant, rrie is the mass of an electron 
and c is the speed of light. Note that the origin of the co- 
ordinate system is defined to be at the center of the NV^ 
defect, in the vicinity of the vacancy (as depicted in Fig. 

H]). Additionally, note that the term J2i ^^0 ^ ^i)^ 
(where e is the electronic charge) quadratic in the mag- 
netic ficld^ has been neglected in the above definition 
since it does not induce a relative shift of the fine struc- 
ture levels or couple the electronic states of the ground 
triplet at first-order in the spin coupling coefficients. The 
second term in ([7]) arises from relativistic corrections to 
the non-rclativistic first tcrm=?- and can be written in the 
more explicit form 
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where the orbital operator G is the matrix 
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The tensor components of the orbital operators con- 
tained in d?]) are 
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Hence, using the general matrix representations of table 
im the matrix representation of the magnetic interaction 
([7]) with the ground triplet can be obtained correct to 
first-order in the spin coupling coefficients 
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where 



91- = .9e + (S2,6 /^)'o,-E + lj^9l,Al + 52,Ai , 

9\\ =5e+5l,Ai 



(11) 



and la,E = -i{ai\\lE\\e)/^/2h, gi^A^ = 

(e||Gi,Ai||e)/2TOeC^ and ^2,^1 = (e||G2,Ai||e)/2meC^. 
Given that the relativistic term ([S]) is much smaller 
than the non-relativistic term, only contributions of the 
relativistic term that are zero-order in the spin coupling 
coefficients have been included. 

Comparing (jlO|) with the matrix representations of 
the total spin operators in table |lll it can be seen that 
the magnetic interaction can be written in the spin- 
Hamiltonian form Vmag = ^S -g ■ B, where the effective 
g-factor tensor g is defined as 



.9± 
g= \ g± 

.gii 



(12) 



The parameters g± and gy may then be identified as the 
non-axial and axial effective g-factor components respec- 
tively. The g-factor components have been measured by 
several ESR studiea^>^i^ and the observed values are 
contained in table IIIII The measurements conclusively 
show that both gy and g± are shifted by +6 ± 1 x 10""*^ 
from the free electron g-factor ge = 2.0023 and one 



TABLE II: Matrix representations of the components of the total spin operator 5* — ^^ Si and tlie orbital tensor operators 
Op,q ~ X^i Op^q{fi) of different symmetry (p, q) correct to first-order in the spin couphng coefficients in the basis of the corrected 
spin-orbit states {'&i°_Ai , ^2°'e,x, ^2°E,y} of the ground triplet. The orbital parameters in terms of reduced matrix elements of 
the center's molecular orbitals are Oa,Ai ~ 2({a;i||C'Ail|ai) + (e||OAil|e)) and Oa,B = (ai||OAi ||e)). 
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study^ observed a small anisotropy of gj^—^ll = 2x10^*. 
From (fTTj) it is clear that there is no orbital magnetic mo- 
ment contribution to gu at first-order in the spin coupling 
coefficients and due to the expected orders of magnitudes 
of the coefficients, any second-order contribution would 
be much smaller than the observed shifts of g|| and g±. 
Thus, only the rclativistic term ([5]) which shifts gn from 
ge by the addition of gi.Ai and also shifts g± by the ad- 
dition of hgi.Ai + 52, Ai I can explain the observed shifts. 
The different contributions of the rclativistic term to g± 
and ^11 can also explain the small anisotropy, however if 
the orbital magnetic moment contribution to g± is large 
enough, it may also contribute at the same order as the 
rclativistic term. The leading order term of the orbital 
magnetic moment contribution to g± is 
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(13) 



due to the expectation that Aj^ ^ vDi-E.2- Through the 
observation of the magnetic circular dichroism (MCD) 
of the optical ZPLf2^ the orbital magnetic moment pa- 
rameter will be of the order la,E ^ 10^^. Consequently, 
the orbital magnetic moment will only contribute signif- 
icantly to the g-factor anisotropy if 2-^^ ^ 10^^, which 

given Ee;i ~ 2.180 eVr^ implies that Aj^ would have to 
be of the order of 1 meV ~ 1 THz, a dramatic differ- 
ence from the axial spin-orbit parameter Ay = 5.3 GHz. 
Such a dramatic difference is improbable given the NV^ 
center's small departure from the higher Td symmetry, 
in which the axial and non-axial spin-orbit parameters 
are equal. Therefore, it appears unlikely that the orbital 
magnetic moment contributes significantly to g. 



B. Interactions with electric fields 

Defining E to be the applied electric field that is as- 
sumed to be approximately constant over the dimensions 



of the NV center, the interaction of the center's elec- 
trons with the electric field is described by the potential^ 



Ki^^rf, 



E 



(14) 



where d ~ efis the electron electric dipole moment. The 
tensor components of d are simply d = dE.x'X + dE.yV + 
d^jZ and the ground triplet matrix representation of Vei 
correct to first-order in the spin coupling coefficients is 
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Vc\ ^ \ d'^E^ da.AiEz + dj^E.j. —d±Ey 
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where da,Ai = 2{{ai\\dAi\\ai) + (e||d^J|e)), dj_ = 
-'2s2,5da,E; d'^ ^ (s2,6 + ^)da.E and da,E = 
7f(ailMB||e). 

The leading order terms of the two non-axial dipole 
parameters are 
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which due to the potentially similar orders of magnitude 
of the numerators in each, suggests that the parame- 
ters arc potentially of the same magnitude. However, 
since d!^ couples electronic states separated in energy by 
Dgsi it can be ignored for static electric fields that satisfy 

d'j_Ej_ <C Dgs, where Ej_ ~ \ E'l + E^ is the non- axial 

electric field strength. This conclusion is in agreement 
with obscrvation,i2i^ where small linear Stark splittings 
of the 771s = ±1 fine structure levels have been shown to 
agree with d^ « and dj_/h = 17±3 Hz cm/VM 

Since correct to first-order in the spin coupling coef- 
ficients, the term da^AiEz is common to each of the di- 
agonal matrix elements, it appears that an axial electric 



TABLE III: Experimental measurements and ab initio calculations of the effective g-factor and hyperfine parameters of the 
NV^ center tabulated by reference. 



Ref. 


9± 


9\\ 


Aj_/h (MHz) 


A^Jh (MHz) 


P/h (MHz) 


Loubseiiii (Exp) 


2.0028(3) 


2.0028(3) 


- 


(±)2.32±0.01 


- 


HeS (Exp) 


2.0028(3) 


2.0028(3) 


(+)2.10±0.10 


(+)2.30±0.02 


(-)5.04±0.05 


Felton^ (Exp) 


2.0031(2) 


2.0029(2) 


(-)2.70(7) 


(-)2.14(7) 


(-)5.01(6) 


Steineii^ (Exp) 


- 


- 


- 


-2.166 ±0.01 


-4.945 ± 0.01 


Smeltzei-i^ (Exp) 


- 


- 


- 


-2.162(2) 


-4.945(5) 


Gali^ {Ab initio) 


- 


- 


(-)1.7 


(-)1.7 


- 



field does not induce a relative shift of the fine structure 
levels of the ground triplet. However, this conclusion is 
in conflict with observation f2£ where a very small linear 
shift of the zero-field splitting between the rrig — and 
rris = ±1 fine structure levels was observed and could be 
described only by a difference of dy/ft, = 0.35 ± 0.02 Hz 
cni/V in the axial dipole parameters of the rus = and 
771s = ±1 spin-orbit states. Such a difference occurs at 
second-order in the spin coupling coefficients, where the 
matrix representation of the interaction of the ground 
triplet with an axial electric field becomes 



d-AiEz = 



where rfii = 
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(17) 
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S2g)d6,Ai and dh,Ai = 
(elld^Jle) — (ailld^Jlai). Therefore, the effective ma- 
trix representation of the interaction of the ground triplet 
with a static electric field satisfying d'j_E± <^ Dgs is 



Vel 




(18) 



This effective representation can be expressed in the spin- 
Hamiltonian form Vel = ^d\\E^Sl - ^dA_E^{Sl - Sl) -f 
-^d±Ey{SxSy + SySx) uscd to describe the linear Stark 
effect present in C^y symmetric systems in ESR.— 

The dipole reduced matrix element da,B contained in 
d± is also responsible for the center's optical transitions^ 
and an estimate of its magnitude can be obtained from 
the center's observed radiative lifetime Tr w 13 ns^ 
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where (Eq) = J^ F{Eo)EQdEo is the expectation value 
of the cube of the optical emission energy given the nor- 
malized vibrational sideband distribution F{Eo)r^ and 
nu = 2.418 is the refractive index of diamond. As the 
vibrational sideband of the center's optical emission ex- 
tends from approximately 1.4 eV to the ZPL at 1.945 
eV, the estimate of the dipole reduced matrix element is 



bounded by 3.65 < da,E/h < 5.98 MHz cm/V (compare 
with, for example, the 5.41 MHz cm/V dipole moment of 
the 5s(25'i/2) ^^ 5p(2P3/2) transition of ^^Rb)^ 

The observed value of d±/h = 17 ± 3 Hz cm/V, ^ 
the estimated range of da,E, the approximate expres- 
sion for d±^ ss '2-^^^da,E, and r/ = 0.053,— imply that 
■^^~ ^ 10^"* as expected, thereby supporting the as- 
sertion made in the previous subsection that the orbital 
magnetic moment does not contribute significantly to g. 
The dipole reduced matrix element dfc,Ai contained in dn 
also contributes to the shift of the center's optical ZPL 
in the presence of an axial electric field<^ However, since 
the optical transition involves a change in MO config- 
uration and, thus a change in the nuclear equilibrium 
coordinates^S both the axial electric and nuclear dipole 
moments contribute to the shift of the center's optical 
ZPL. Consequently, without knowledge of the nuclear 
dipole moment, it is not possible to estimate d^.Ai given 
just measurements of the shift. 



C. Interactions w^ith strain fields 

The interaction of the center's electrons with a crystal 
strain field can be approximately described by perform- 
ing a Taylor series expansion of the electronic Hamil- 
tonian He in terms of the displacements of the nuclear 
coordinates R from their ground state equilibrium coor- 
dinates Rq induced by the strain field and retaining only 
the linear terms of the expansion. Defining Qu,p,q to be 
the u*^ normal nuclear displacement coordinate of the 
crystal that transforms as the row q of the irreducible 
representation p of the C^y group, the strain potential 



^- = EE 



dVNe{n,R) 



I u,p,q 



dQu. 



p,q 



u 



p,g 



(20) 



where ^u,p,q is the crystal strain along the displace- 
ment coordinate Qu,p,q- Given that by definition 
dyNe(fi,R)/dQu,p,q\ii is an orbital tensor operator of 
symmetry [p, q) , the ground triplet matrix representa- 
tion of Vstr is analogous to that of Vei- Treating the 
non-axial terms correct to first-order and the axial terms 
correct to second-order in the spin coupling coefficients. 



the matrix representation of Vstr is 

C C C 

v.tr =\ C C + C- + C. ^ -C, 

Cj/ "^y ^z ~^ ^z ~ ^x 



(21) 



where 

c 

Cz 

CI 

a 



u 

(*2,5 + *2,6 + ■52,9) 2^ Cu,b,Ai^u,Ai 
u 

(S23 H ■j=) 2_^ (uM,E^u,E,k 

~2S2,5 ^_^ Cu,a,E^u,E,k 

dVNe ,, , .. ,, ay, 



= 2 a 



(el 



dQu,Ai 

dVNe 



ai) +2(e 



TVe 



5Q«,Ai 



(«l|| 



i9Qu,Ai 

dVNe 



|e) 



(fli 






dQu,Ai 



m) 



Ro 



V2' " dQu. 



E 



(22) 



-Ro 



and k = X, y. Note that torsional strain components that 
have displacement coordinates that transform as A2 have 
been ignored. 

Similar to Vei, the diagonal contributions C^ do not 
shift the fine structure levels of the ground triplet with 
respect to each other and thus can be effectively ignored. 
Likewise, for non-axial strains that satisfy ^^ ^ Dgs 



(where C,'^ = J C,'^ + Cy'^); Ci and C,y will have negligible 
effect on the fine structure levels and state coupling and 
thus can also be ignored. Defining the effective strain 
field a = aa:X + ayy + azZ = Cx/d±x + Cy/d±y + Cz/d\\z, 
the effective matrix representation of Vstr becomes 





K;tr = I d||(T^ +d±ax 

—d±ay 




(23) 



and it is clear that the strain field can be treated as 
an additional effective local electric field <t at the cen- 
ter. Hence, by defining the total effective electric field 
II = E + (T, the interaction of the center with both 
strain and electric fields can be expressed in the spin- 



Hamiltonian form Vei + Vstr = ■p-rfyn^S'^ - 

y) ~^ J^dj_Ily[DxSy + OyOx)- 



w 



d±UxiSl^ 



D. The complete spin-Hamiltonian 

Including the descriptions of the spin-spin zero-field 
splitting and interactions with magnetic, electric and 
strain fields obtained in the previous subsections, the 



complete electronic spin-Hamiltonian of the ground state 
spin becomes 



H„ 



~T^d±iix[,^x ^ ^y) ^ '7^d±iLy\yOxOy + Oybx) 

(24) 

which in the spin basis {SA2i-'SE.y, Se,x} associated 
with the ground triplet spin-orbit states, has the matrix 
representation 



H, 



gs 




-iBy 



iBy 



-iBx 



V + Ex —iBz — Ey 



(25) 



iBx iBz — £y T> — £x 



where V = Dgs + d||H2, B^ = fj,Bg\\Bz, £k = d±Ilk, Bk = 
fJ-Bg±Bk, and k=x,y. This final matrix representation 
provides the simplest description of the dependence of 
the ground state spin on the six independent electric- 
magnetic-strain field parameters (2?,fa;,fy,K) and will 
be used to obtain the spin solution in part II of this 



paper series 
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III. NUCLEAR HYPERFINE STRUCTURE 

The nuclear hyperfine interaction Vhf = Vmhf + Vehf 
between the center's electrons and the nuclei of the crys- 
tal lattice has magnetic Vmhf and electric T^^/ compo- 
nents. The magnetic hyperfine component accounts for 
the interactions of the electronic spin and orbital mag- 
netic moment with the nuclear spins of the lattice. Since 
it was found in the previous section that the ground 
triplet has no orbital magnetic moment at zero-order in 
the spin coupling coefficients, the interaction between the 
orbital magnetic moment and the nuclear spins maybe ig- 
nored. The electric hyperfine component approximately 
accounts for the finite size of the nuclei through the 
interaction of the center's electrons with the electric 
quadrupole moments of the finite charge distributions 
of the nuclei. In an isotopically pure ^^C crystal, only 
the ^*N nucleus of the center will have a non-zero nu- 
clear spin (/ = 1) and a non-zero electric quadrupole 
moment jS^ Consequently, in an isotopically pure crystal 
the magnetic and electric hyperfine interactions of the 
NV~ center are described by the potentials^ 



Vr 



mhf 



Cmhf Y, ( 47r<5(f,w) - -7-^ ) s, ■ I 



riN\ 



3(si ■fiN){riN ■ I) 



\riN\ 



V 



ehf 



2Zn -^, ^ 



d'VM) 



N ' — ', \ ^-^ dRN.adRN,b 



Ro 



d^V 



dRN,adRN,b 



^a,a^a,b 



(26) 



where C^hf = fJ-BfJ-NgegN^j^, fJ-N is the nuclear mag- 
neton, 5JV = 0.40356 is the ^''N nuclear g-factor,— / is 
the spin operator of the ^''N nucleus, fiN = Rn — ^i ^ 
XiNX + yiNV + ZiNZ, Rn = Rn,xX + RN,yy + Rn,zZ is the 
position of the ^^N nucleus, Ve is the Coulomb interaction 
potential of the center's electrons with the ^^N nucleus, 
Vi is the Coulomb interaction potential of the lattice elec- 
trons and "'^^C nuclei with the ^^N nucleus, Zn = 7 is the 
relative charge of the ^''N nucleus, and Ua^a and Ua^b are 
the components of the displacement of the a*'' ^''N pro- 
ton from Rn in the a,b = x,y^ z directions. 

The magnetic hyperfine interaction can be written as a 
sum of terms containing rank two orbital tensor operators 

V,nhf =J2i^i' i^AiirrN) + AE,x{,rrN) + ^B,y(nw)] ' I- 

By applying the matrix representations of table |lll it is 
clear that only the term containing Aai (fjAr) contributes 
at zero-order in the spin coupling coefficients. Given this 
result, Vrnhf can be written in the spin-Hamiltonian form 
Vmhf ^ S ■ A- 1, where 



A = 



A_L 
A_L 

All 



A\\ = /ai + 2aAi is the axial magnetic hyperfine 
parameter, Aj_ = /ai — ciai is the non-axial mag- 
netic hyperfine parameter, Jai = C'm/i/47r(e||(5(fijv)||e) 
is the Fermi contact contribution, and uai — 

kCrr^^hfieW-r-^ (ifep' ~ 1) ll*^) '^ *^® dipolar contribu- 



■^mhf\ 



tion. Measured values of the magnetic hyperfine param- 
eters are contained in table Hill and, although they differ 
in sign and magnitude, conclusively show that both the 
Fermi contact and dipolar contributions must be non- 
zero. For example, using the values obtained in Ref. |44| . 
f^Jh = -2.51 MHz and uaJH = 187 kHz. 

The expression for the Fermi contact contribution can 
be simplified further to /ai = C'„i/i/47r|ea;(i?Ar)p = 
Cmhf4:TT\ey{RN)\'^- As the e MOs transform as the E 
irreducible representation of the Csv group, they are by 
definition zero at any point along the axial symmetry 
axis of the center, and since the equilibrium position of 
the ^^N nucleus is on the axial symmetry axis, the Fermi 
contact contribution vanishes if the ^*N is fixed at its 
equilibrium position. To account for the non-zero Fermi 
contact contribution in the molecular model, the vibra- 
tional wavefunction xn{Rn) of the ^^N must be consid- 
ered, in which case the expression for the Fermi contact 
contribution becomes 



Jvib — ty, 



nihf 



Att 



\eARN)r\XN{RN)rd-'RN (27) 



A similar vibrationally corrected expression for the dipo- 
lar contribution can also be defined. 

The magnetic polarization of the ^"^N core electrons 
associated with the nis = ±1 states of the center will 
provide a negative Fermi contact contribution /eor d^^'^^ 



in addition to the positive contribution arising from the 
vibrational motion of the nucleus. However, the descrip- 
tion of fcore is beyoud the molecular model in its current 
formulation, as only the interactions of the bound va- 
lence electrons are considered in the model. An ab initio 
study^ has yielded fcore ~ —1.7 MHz, but no ab initio 
calculation of the vibrationally corrected /vib has been 
conducted to date. Although fmi, is likely to be much 
smaller than fcore, the calculation of /vib will provide a 
more complete ab initio model of the magnetic hyper- 
fine parameters. Likewise, an extension of the molecular 
model to describe fcore will also provide further insight 
into the interactions between the nucleus and the bound 
electrons. 

The electronic component of the electric hyperfine in- 
teraction Vehf can be written as a sum of products of 
orbital electronic and nuclear tensor operators. Anal- 
ogous to the interaction of the ground state spin with 
electric and strain fields, at zero-order in the spin cou- 
pling coefficients only the terms of Vehf that contain Ai 
orbital electronic tensor operators will have non-zero ma- 
trix elements, and these matrix elements will be diagonal 
and identical for each spin state. The symmetry of the 
lattice ensures that only similar terms from the lattice 
electron and internuclear contribution are also non-zero. 
Consequently, the electric hyperfine interaction within 
the ground state spin reduces to 



V, 



ehf 



4Z 



N 



E3^ 



(28) 



2{e\\d'^Vc/dR 



'N,z\Rq\ 



ziflo 



\ai) 



d^Vi/dRlJ^^ 



IS pro- 



where q^ 

2{a^\\d^Vc/dRl^ 

portional to the axial gradient of the net axial electric 
field at the -'^'^N nucleus generated by all of the electrons 
and other nuclei. 

As per standard practice, the nuclear quadrupole op- 
erator X]q 3"a,z ~ l^ap can be replaced by a nuclear 
spin operator Qzl^/f^ of the same symmetry through 
the definition of the proportionality constant Qz = 
(1' 1| Y.a ^^a.z - |"aP|l, 1) usiug thc 4 dgcustate |/ = 
I, nil ~ l)'^ Qz therefore quantifies the difference in 
the axial anisotropy of the nuclear quadrupole moment 
between the mj = ±1 and mi = states. The final 
form of the effective electric hyperfine interaction is then 
Vehf = Plz/f^^j where P = QzQz/^Zn is the nuclear 
quadrupole parameter contained in table IIIIl 

Combining the expressions obtained for the magnetic 
and electric hyperfine interactions, the zero-field spin- 
Hamiltonian of the ground state spin including nuclear 
hyperfine interaction is 



Hhf = ^[D,,Sl + A^\SJz 



+ A^{S^h+SyIy)+Pll] 



(29) 



The approximate solutions of Hhf can be obtained by 
constructing thc hyperfine states of the ground triplet 



in a similar manner to the earlier construction of the 
electronic spin-orbit states by defining linear combina- 
tions of products of electronic spin-orbit and nuclear spin 
states that have definite Cs^ symmetry. Given that the 
symmetrised nuclear spin states of the ^'^N nucleus in 
terms of the Iz eigenstates {|/, m/)} are IA2 ~ II7O), 
Ie,. = ^(|l,l) + |l,-l))and/B., - ^(|1,1)-|1,-1)), 
the symmetrised hyperfinc states '^n.j,k arc 



*3,B,y = -^{^tsjE^y - ^2,E..y 



*4.Ai-^($2°B,x^£,^ + * 



Ie^x) 

2,E,ylE,y) 



*5.A, = -^{^EjE^y 



*, 



Q,E,x 



-$ 



Ja2, *6 



KE,ylE,. 



2,E,y^A2^ '*'6,E,y — ^2,E,x^A2 



(30) 



^l,E,x 
^2.A2 -- 



'^T^aJe,. 

nyjA. 



*l,£.y = ^T,A^lE,y 



*3,£,. = -j^m:E,xiE,x - n:E,yiE, 



The matrix representation of the zero-field Hamilto- 
nian in the basis of hyperfine states {^2,^2: ^q.e.xi 

"^Q.E.y, *l,_E,i^, *l,£,y, *3,£,x, *3,B,y, ^4,Ai , ^S.As} is 
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(31) 



r 



where h± = P ± A||. Note that the above basis of hy- 
perfinc states has been grouped into states of the same 
nuclear spin projection. The energies E^^ correct to 
first-order in nuclear hyperfinc interactions can be eas- 



ily inferred: E^ 



hf 



e: 



hf 



E, 



hf 



D 



gs 



P,E'2 

p 



hf 



0, ^,f 



D 



An, E, 



hf 



gs 



P + An 



Dgs] and the 



corresponding hyperfine structure is depicted in Fig. |3l 
Since A± is observed to satisfy A± ^ Dgs ^^^ the non- 



axial magnetic hyperfine parameter that couples hyper- 
fine states of different nuclear spin projection will have a 
negligible effect on the zero-field hyperfine structure and 
state couplings. 

The matrix representation of the electronic interaction 
with electric, magnetic and strain fields in the hyperfine 
basis is 



KiaK + V^l + y,tr 
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(32) 



The matrix representation demonstrates that if the much smaller interaction of the nuclear spin with the fields is 
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FIG. 3: (color online) (a) The hyperfine structure of the ground triplet. The level splittings are indicated by dashed arrows 
and the observed values for the hyperfine parameters Ay and P are contained in table Hill The ordering of the hyperfine levels 
is based upon the parameters measured in Ref. |4J- The allowed magnetic transitions between levels of the same nuclear spin 
projection are indicated by solid arrows and labeled by their transition energy AE^f^. (b) The splittings of the rus = ±1 
hyperfine levels in the presence of non-axial electric-strain fields of magnitude 11 _l (upper) and axial magnetic fields Bz (lower). 
mi — hyperfine levels are colored blue and mi = ±1 hyperfine levels are colored red and green. The hyperfine parameters of 
Ref. |4J have been used to calculate the splittings. 



ignored, the fields do not couple states of different nu- 
clear spin projection. Furthermore, by comparing the 
above matrix representation with that of Hgs, the rep- 
resentation also demonstrates that the hyperfine states 
of nuclear spin projection mj = (upper 3x3 diago- 
nal block) interact with the fields in a manner similar 
to the electronic spin-orbit states discussed in the pre- 
vious section, whereas the to/ = ±1 hyperfine states 
(lower 6x6 diagonal block) interact differently in the 
weak field limit, where the fields induce shifts compara- 
ble to the hyperfine splittings. For example, as depicted 
in Fig. [31 the (rUs = ±1, to/ = 0) states split linearly in 
the presence of a non-axial electric-strain field, whereas 
the non-degenerate (tos = ±1, mj = ±1) states repel 
quadratically in the presence of non-axial electric-strain 
fields that satisfy £± < 2A|[ . Note that in the large field 
limit, where the fields induce shifts much larger than the 
hyperfine splittings, both sets of hyperfine states behave 
approximately analogous to the spin-orbit states. 

The allowed magnetic transitions are depicted in Fig. 
13] and indicate that in the absence of static fields, there 



will exist three lines in the hyperfine spectra with ener- 



gies Ae!;{ 



D„ 



- A|| , AE'1{ 



ae!^!i 



Dgs — j4|| and 



AEq 2 = Dgs in agreement with observation!^ The cen- 
tral hyperfine line therefore corresponds to transitions 
between mj = states and the lower and higher en- 
ergy lines correspond to transitions between mi = ±1 
states. Consequently, due to the different interactions of 
the To/ = and to/ = ±1 states in the weak static field 
limit, the central hyperfine line will depend differently 
on the static fields compared to the lower and higher en- 
ergy lines. These differences in the dependence of the 
hyperfine lines in the weak field limit was used in the re- 
cent electric field sensing demonstration)^ where the mag- 
netic field was precisely aligned in the non-axial direction 
(Bz = 0) in the presence of a non-axial electric-strain 
field by observing the splitting of the mi = ±1 hyper- 
fine lines whilst the measurement of the electric field was 
conducted by observing the linear splitting of the central 
TO./ = hyperfine line. Hence, it is clear that the hyper- 
fine structure of the ground state spin and its more com- 
plicated interactions with electric, magnetic and strain 
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fields, is an important consideration for applications of 
the spin that operate in the weak field limit. 



IV. CONCLUSION 

In this article, the theory of the ground state spin has 
for the first time been fully developed using the molecu- 
lar model of the center in order to provide detailed expla- 
nations for the spin's fine and hyperfine structures and 
its interactions with electric, magnetic and strain fields. 
Given these explanations, an effective spin-Hamiltonian 
that describes the electronic states in the high field limit 
and the to/ = subset of hyperfine states in the low field 
limit was derived. The explanations also allowed the cor- 
relation of the properties of the ground state spin with 
the other properties of the center and provided explicit 
expressions for the key parameters of spin in terms of 



the center's MOs. Hence, this work has identified the 
critical parameters that need to be pursued by future ex- 
perimental and ab initio studies. Furthermore, this work 
has also provided the essential theoretical understanding 
of this remarkable spin that can be used to model the 
spin in its ground-breaking quantum metrology and QIP 
applications. 
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